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Abstract 

The nature of current sheet formation in the vicinity of three-dimensional (3D) magnetic null 
points is investigated. The particular focus is upon the effect of the compressibility of the plasma 
on the qualitative and quantitative properties of the current sheet. An initially potential 3D null 
is subjected to shearing perturbations, as in a previous paper [Pontin et al, Phys. Plasmas, in 
press (2007)]. It is found that as the incompressible limit is approached, the collapse of the null 
point is suppressed, and an approximately planar current sheet aligned to the fan plane is present 
instead. This is the case regardless of whether the spine or fan of the null is sheared. Both the peak 
current and peak reconnection rate are reduced. The results have a bearing on previous analytical 
solutions for steady-state reconnection in incompressible plasmas, implying that fan current sheet 
solutions are dynamically accessible, while spine current sheet solutions are not. 
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I. INTRODUCTION 



In astrophysical plasmas, such as the solar corona, the three-dimensional (3D) magnetic 
field topology is often highly complex. In such complex 3D magnetic fields, where traditional 
two-dimensional (2D) X-point magnetic reconnection models may no longer be applicable, 
determining the sites at which dynamic phenomena and energy release may occur is a crucial 
and non-trivial problem. Due to the typically very high Lundquist number, such events occur 
only at locations where intense currents (singular under an ideal MHD evolution) may form. 



One such site is a 3D magnetic null point (e.g. Refs. 
formation at such 3D nulls is investigated here. 
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3D null points are predicted to be present in abundance in the solar corona (e.g. Refs 
Q]). Furthermore, there is observational evidence that reconnection at a 3D null may be 
important in some solar flares^, as well as in eruptive phenomena in active regions^. In 
addition, the first in situ observation^ of reconnection occurring at a 3D null point in the 
Earth's magnetotail has recently been made by the Cluster spacecraft. Moreover, current 
growth at 3D nulls has been observed in the laboratory^. 

The magnetic field topology and geometry in the vicinity of such a null can be described 
by the two sets of field lines which asymptotically approach, or recede from, the null. A pair 
of field lines approach (recede from) the null from opposite directions, defining the 'spine' 
(or 7-line) of the null. In addition, an infinite family of field lines recede from (approach) 
the null in a surface known as the fan (or S-) plane (see Refs. 

To this point, many studies of the MHD behaviour of 3D nulls have been kinematic, see 
e.g. Refs. jQQllsj]- However, a few solutions to the full set of MHD equations do exist 
for reconnection at current sheets located at 3D nulls, in incompressible plasmas. These 
incompressible solutions are based upon the technique first proposed by Craig & Henton^ 
for the 2D reconnection problem. The solutions describe steady-state current sheets aligned 
to the fan and spine of the nult^. Time-dependent solutions for the fan current sheets also 
exisir^. 

In a previous paper — Ref. hereafter referred to as paper I — we investigated the be- 
haviour of 3D null points which are subjected to shearing boundary motions, and found that 
current sheets formed at the null. In this paper we consider the effect of moving from the 
compressible towards the incompressible limit, which is found to have a profound effect on 
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both the quantitative and qualitative properties of the current sheet. This is highly rele- 
vant when it comes to comparing the observed current sheet formation with the analytical 
models, which must by necessity invoke various simplifications. Typically, the plasma in 
the solar atmosphere or Earth's magnetosphere is compressible. Thus it is of great interest 
to understand the relationship between this regime and the incompressible approximation, 
upon which much of the previous theory has been based. 

The remainder of the paper is set out as follows. In Sec. [Ill we briefly review the pre- 
vious results of paper I. In Sec. II III we describe simulations in which we move towards the 
incompressible limit, and in Sec. IIVI we discuss the relation of our results to analytical 
incompressible solutions, and the implications for their dynamic accessibility. In Sec. IVl we 
consider the case where we drive across the fan instead of across the spine of the null, and 
finally in Sec. [VI] we present a summary. 



II. BEHAVIOUR IN A COMPRESSIBLE PLASMA 



In paper I, we discussed the evolution of the magnetic field in the vicinity of a generic 
3D magnetic null. We demonstrated by means of a kinematic solution that an evolution of 
the null which acts to change the angle between the spine and fan (such that the ratios of 
the null eigenvalues change in time) is prohibited in ideal MHD. We then went on to present 
the results of numerical simulations, which demonstrated the formation of strong current 
concentrations at the null in response to boundary perturbations. Simulation runs based on 
the same numerical code are presented below (for further details on the numerical scheme, 



see Refs. 
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20|). 



At t = the magnetic field in the domain is given by B = B (—2x, y, z), which defines a 
3D null whose spine lies along the x-axis, and whose fan is in the x = plane. J = 0, and so 
taking the density (p) and internal energy (e) of the plasma to be uniform at t = we begin 
with an equilibrium [we take p = 1, e = ^7/(7 — 1), where f3 is a constant which determines 
the plasma-/? (which is of course spatially dependent) and 7 is the ratio of specific heats]. 
All of the domain boundaries are line-tied, and are located at [x,y,z] = [±Xj, ±YJ, ±2Tj]. 
The configuration is then perturbed by imposing a plasma flow on the x-boundaries, while 
the y- and z-boundaries are placed sufficiently far away that there is insufficient time for 
information to propagate to them and back to the null before the simulations are halted. 
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A single time unit in the simulations is equivalent to the Alfven travel time across a unit 
length in a plasma of density p = 1 and uniform magnetic field of modulus 1. The resistivity 
is taken to be uniform, with its value being based upon the dimensions of the domain. Note 
that at t = 0, B is scale-free as it is linear, and thus, the actual value of t] is somewhat 
arbitrary until we fix a physical length scale to associate with the size of our domain. 

The boundary driving takes the form on each boundary of two distorted vortices of 
opposite sense, which combine to provide the desired effect of advecting the spine in the y 
direction, in opposite senses on opposite boundaries (x = ±Xi) [see paper I, Eq. (19) and 
Fig. 2(b)]. In the majority of the runs described, the driving profile is transient, with its 
time dependence defined by 



Vq, t constant, so that the driving switches on at t = and off again at t — 2r. The result is 
that a current concentration forms at the null, which is expected to be singular in the ideal 
limitA During the early evolution, a stagnation flow, accelerated by the Lorentz force (but 
opposed by the plasma pressure) acts to close up the spine and fan towards one another 
locally at the null. The initial null is unstable to such a collapse of the spine and fan in 
any plane containing the spine, with the z = plane being selected by the orientation of 
the boundary driving. It is precisely this collapse that was shown in the kinematic solution 
earlier in the paper to be prohibited under ideal MHD. Thus it must be facilitated by 
non-ideal processes. 

Due to this local collapse, a current sheet forms which typically spans the collapsed spine 
and fan, with a tendency to spread along the fan surface (especially for weaker driving). 
Accompanying the current growth is the development of a component of E parallel to B 
(E»), signifying a breakdown of ideal behaviour, and magnetic reconnection. The integral 
of this quantity along the magnetic field line in the fan perpendicular to the shear plane can 
be shown to give a physically meaningful measure of the reconnection rate — giving the rate 
of flux transfer across the fan (separatrix) surface, see Ref. [l^]. 

An examination of the quantitative properties of the current sheet showed that the peak 
current, peak reconnection rate, and sheet dimensions all scale linearly with the modulus 
of the driving velocity. In addition, under continual boundary shearing, the current sheet 
appears to grow in size and modulus indefinitely (rather than being controlled by any self- 




< t < 2r, 



(1) 
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regulating mechanism). This type of behaviour is also observed in 2D 'forced' or 'driven' 
reconnection simulations. The nature of the current sheet seems to be controlled at all 
times by the degree of boundary displacement of the spine and fan (and how quickly this 
displacement is attained), and so it may share some properties (at a given time) with the 



2D 'non-uniform reconnection' regimes (see Ref. [2l| for a review). Care should be taken, 
however, in drawing parallels with either of these 2D models, since each involves an inflow 
of plasma through the boundaries. By contrast, in our simulations the driving velocity is 
imposed parallel to the boundaries. 

In paper I we considered the case of a monatomic ideal gas, that is we took the ratio of 
specific heats, 7 = 5/3. It is straightforward to see that the incompressible limit may be 
reached formally by letting 7 — > 00. Taking the time-derivative of the polytropic equation 
of state, p/p 1 = const, and substituting for dp/dt using the continuity equation gives 

1 dp 

V ■ v = — . 

7P at 

III. TOWARDS INCOMPRESSIBLE LIMIT 



We repeat here the simulations described in paper I, with increased values of 7. This is 
somewhat problematic numerically (due to the increased wave speeds in the system), but in 
fact it turns out that even for moderately large values of 7, the differences are striking. 



A. Qualitative differences 

The parameters chosen for the simulation runs closely follow those taken in paper I, and 
are as follows. We take B = 1, the driving strength w = 0.01, r = 1.8, = 80 (boundary 
driving localisation), (3 = 0.05, 77 = 5 x 10 -4 , and the numerical domain has dimensions 
Xi = 0.5, Yi = Z[ = 3. The resolution of the simulations is 128 3 , on a non-uniform mesh 
with smallest grid spacing near the null to achieve higher resolution there; 5x ~ 0.0035 and 
6y,6z ~ 0.020. 

As the driving begins (t = 0), a disturbance propagates along the spine (and nearby 
field lines), and focuses at the null. For 7 = 5/3, the null point 'collapses' with the spine 
and fan closing up towards one another. A strongly focused current sheet spans the spine 
and fan. However, for larger values of 7, the current concentration distributes itself along 
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FIG. 1: Isosurfaces of |J| at 50% of maximum, at the time of its temporal peak, for (a) 7 = 5/3, 
(b) 7 = 10/3 and (c) 7 = 10. 






-0.4 < x >0.4 -0.4 < x >0.4 -0.4 < x > 0.4 

FIG. 2: (Colour online) Current density |J| in the z = plane, at the time of its temporal peak, 
for (a) 7 = 5/3, (b) 7 = 10/3 and (c) 7 = 10. 
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the fan surface, becoming more weakly focused at the null for increasing 7 (see Figs. HJ [2]). 
Furthermore, the fan surface remains increasingly planar at larger 7 (see Fig. [2]), and also the 
spine and fan do not collapse towards each other to the same extent. This is demonstrated 
in Fig. H](f), where the minimum angle between the spine and fan (9 m i n ) is plotted for runs 
with various values of 7. We observe that even for 7 = 10, although the current sheet is 
approximately planar (at x = 0), the minimum angle between the spine and fan is still 
significantly less than ir/2. This is because the spine is still driven towards the fan by the 
boundary driving (most of the stress from which is taken up in the weak field region around 
the null itself), even though the fan remains approximately in the x = plane rather than 
collapsing sympathetically towards the spine. 

It is worth noting that the above described behaviour also depends on other parameters 
in the simulation. For example, how effectively the null collapses is also dependent on the 
driving speed, with greater collapse and stronger focusing of the current sheet for larger 
vq (see paper I). Therefore larger values of 7 are likely to be required in order to render 
the fan approximately planar for larger Vq, and also for larger r (longer driving time). The 
plasma-/? is also a crucial parameter, since we find that increasing f3 has a very similar effect 
to increasing 7. It is natural to expect this on physical grounds, since increasing either 
parameter has the effect of increasing the sound speed, and reducing the effect of magnetic 
forces in plasma compression. Finally, since the null collapse is driven by the Lorentz force, 
a thinner more intense current sheet, which will form for a lower value of rj, will increase 
the degree of collapse. Thus, the extent to which the null collapses and the current focusses 
at the null is dependent on a combination of the driving velocity (i> , r) and the plasma 
parameters 7, (3 and 77. 

An obvious question when examining the above results is whether the planar current 
sheet in the fan plane for large 7 is a result of the symmetry of the configuration, with the 
null at the centre of the domain and the fan plane parallel to the driving boundaries. We 
therefore re-ran the simulations at large 7 with the null point rotated by a finite angle in 
the xy-p\ane (so that the spine and fan were no longer parallel to the boundaries). In this 
case, a planar current sheet still forms in the fan, and thus our results seem general in this 
respect. 

Accompanying the changing current localisation as we move towards the incompressible 
limit is a change in the behaviour of the plasma flow. This again signifies the fact that the 
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FIG. 3: (Colour online) Plasma flow in the z = plane at t ~ 2.5, for (a) 7 = 5/3 and (b) 7 = 10. 
Background shading shows |J|. 

fan of the null remains increasingly planar. For 7 = 5/3, a stagnation flow is typically set 
up, which is accelerated by the Lorentz force (and opposed by the plasma pressure gradient), 
and which closes up the spine and fan. However, for larger 7 this flow is absent, and instead 
v x is approximately zero, and the flow is roughly parallel with the driving boundaries (see 

Fig. ED. 

Finally, it should be noted that all of the above considerations are the same as for the 
case of a 2D X-point. That is, repeating the above simulations but with the magnetic 
field at t = defined by B = B (— x,y,0), we see the same trend. For 7 = 5/3 the X- 
point collapses, forming a current sheet which locally spans the two separatrices (a 'Y-point' 
appearance), but for large 7 the X-point c ollap se is suppressed, and the current spreads 



along the (unsheared) separatrix (as in Ref. [l6j|). 



B. Quantitative differences 

It is not only the qualitative properties of the current sheet which are affected by changing 
the plasma compressibility. Accompanying the spreading of the current sheet along the fan 
for increased 7 is a decrease in the peak current and reconnection rate in the simulation, see 
Fig. H(a, b). The rate of change of each quantity around 7 = 5/3 is much greater than that 
around 7 = 10, implying that even for this moderate value of 7, the behaviour is already a 
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FIG. 4: Scaling with 7 of (a) the peak current (J max ), (b) the peak reconnection rate (J Eu), (c, 
d, e) the current sheet FWHM at time of peak current (L x , L y , L z ) and (f) the minimum angle 
between the spine and fan (6 m in)i where the dashed line indicates the value 6 = ir/2. 

fairly good approximation to the incompressible limit (for all other parameters fixed). The 
change in geometry of the current sheet is evidenced by the variation in the dimensions of 
the region of high |J|, L x , L y and L z (measured at the time of current maximum, by the full- 
width-at-half-maximum (FWHM) in each coordinate direction). L y and L z increase with 7, 
showing how the current spreads along the fan surface as we move towards the incompressible 
limit [Fig. H](d, e)]. On the other hand, L x decreases as 7 increases, demonstrating that the 
null point collapse is inhibited [Fig. H(c)]. Even for 7 = 20, L x essentially measures the 
current sheet 'thickness'. 
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IV. RELATION TO ANALYTICAL SOLUTIONS 



A. Dynamic accessibility 

We now investigate the relation between our simulation results and previous analytical 
solutions for incompressible plasmas. In the steady-state solutions of Craig et al— ^ the 
assumption of incompressibility leads to a symmetry between B and v in the MHD equa- 
tions. Progress is then made by defining a 3D current-free 'background field', upon which 
disturbance fields of low-dimensionality are super-imposed. This necessarily results in cur- 
rent sheets which are also of reduced dimensionality. The solutions are sometimes referred 
to as 'reconnective annihilation—, since they contain current sheets of infinite extent in at 
least one direction, and as a result the plasma advects field lines across either the spine or 
the fan, but they only diffuse towards the other of these (through the current sheet). It 
might be expected that the infinite nature of the current sheets is due to the severe analyt- 
ical restriction of low-dimensionality 'disturbance fields'. However, as we have seen above, 
applying shearing boundary motions to the spine footpoints of the null indeed results in a 
quasi-planar current sheet in the fan plane, albeit only for large 7. 

Of great importance for any steady-state solution is its dynamic accessibility under a 
time-dependent evolution. Investigations into the dynamic accessibility of two-dimensional^ 



solutions have been carried out by various authors (e.g. Refs. 23IJ24I ]). The results of the 
previous section provide strong evidence that in a fully dynamic and fully 3D (yet incom- 
pressible) system, the fan current sheet solutions are indeed dynamically accessible. One 
further question which presents itself here is whether in fact the spine current solutions are 
also dynamically accessible. In the analytical solutions, a tubular spine current results from 
shearing perturbations of the fan plane. This is investigated in Section |V] 



B. Breakdown of analytical solutions 

It appears that in the incompressible limit, fan current solutions are dynamically accessi- 
ble, and (at least qualitatively) provide an accurate snapshot of the dynamical and fully 3D 
behaviour. However, in the case of a compressible plasma this appears not to be the case. 
In order to understand why this is, we must examine the force balance which exists in the 
analytical solutions. 
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The method of the analytical solutions is based upon taking the vector product ('curl') 
of the momentum equation, and solving this in conjunction with the induction equation. 
The pressure can then be calculated a posteriori. However, it has been realised^^^L^ 
this places a limit on the maximum current (or reconnection rate) which can be attained in 
these 'flux-pile- up' solutions, since the current sheet must be maintained by a large pressure 
at infinity. For current values above some limit, the pressure required is unphysically large. 

We can similarly examine the plasma pressure (or pressure gradient) which exists within 
the current sheet itself. In the steady-state fan current solution of Craig et al.—, the magnetic 
and velocity fields are defined by 



A, «, a constant, < k < 1. The pressure is found from the momentum equation, and the 
pressure gradient perpendicular to the fan plane is given by 



Solving the induction equation for Y and Z (see Ref. [30|]) reveals that in the current sheet, 



gradient which scales as a negative power of 77, which becomes extremely large at realistic 
values of rj for astrophysical plasmas. Note though that the strongest pressure restriction 
occurs in the degenerate 2D case (ft = or k = 1). Once the pressure gradient can no 
longer accommodate the huge Lorentz force within the sheet, the null point will begin to 
collapse, and the strict planar nature of the fan plane and current sheet will be lost (note 
that the Lorentz force always points in the direction which further closes the angle between 
spine and fan, while the pressure gradient acts in the opposite sense). With the symmetry 
of the system broken, the analytical solutions can no longer describe the behaviour, and we 
can expect the nature of the current sheet to be significantly altered. A similar argument 
has been made by Ma et al— for the case of disturbances perpendicular to a 2D planar 
X-point — they found that once the strict symmetry of the system was broken, qualitatively 
very different behaviour resulted. 

Leaving the steady-state solutions and examining instead the time-dependent fan current 
sheet solutions^, one arrives at a similar conclusion. In this case, the time-dependent 



B = \P + Y(x)y + Z(x) 



z, v = P + AF(x)y + AZ(x)z, 



P = a (—x, rzy, (1 — k)z) . 




dY/dx ~ rj ( 1+K )/ 2 ; OZ/dx ~ rj ^ K ^ 2 . Thus in the current sheet we require a pressure 
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FIG. 5: Pressure gradient at the time of maximum current in the z = plane, for the run with 
7 = 10, with driving across the spine. 



pressure gradient force in the x-direction in the ideal localisation phase is given by 



for one disturbance component (a~ is a constant which determines the relative strengths 
of the background magnetic and plasma flow fields). This peaks once resistive dissipation 
becomes important and the current density reaches a maximum value, when we have 



The contribution of the other disturbance component is obtained by replacing k by (1 — k) 
in each of the above, < k < 1. Thus the plasma pressure force in the re-direction (or 
symmetry-breaking direction) increases exponentially in time, in order to counteract the 
effect of the increasing Lorentz force. For sufficiently small 77, the pressure force will no longer 
be able to balance the Lorentz force during this localisation process, and the symmetry of 
the configuration will be lost. 

The effect of the pressure gradient within our simulations is shown in Fig. Here, vectors 
of Vp are plotted in the z = plane for 7 = 10 at the time of the peak current. It is clear 
that the pressure gradient force behaves exactly as described — its effect is localised primarily 
within the current sheet (near the x = plane; compare with Fig. [2]), and is directed in such 
a sense as to oppose the collapse of the fan surface and current sheet. 



^1 „ _ e «" (i+nyt 



dx 
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The fact that the geometry of the current sheet which we observe in our compressible 
simulations is very different to that of the analytical solutions is not completely unprece- 
dented. In fact, in laboratory experiments examining the formation of current sheets at 3D 
nulls, Bogdanov et ai— made a similar observation. They too found a current sheet forming 
at a finite angle to the global directions of both the null spine and fan, which had not been 
expected from prior self-similar analytical solutions^. However, it is interesting to observe 
that the incompressible solutions are indeed recovered in the limit of large 7, even though 
we make no assumption regarding the dimensionality of any fields in the solution. 

Note finally that all of the arguments given above carry through to the 2D case. Thus 
our results for the 2D null, when compared with the solution of Ref. 
by similar reasoning. 



16|, can be explained 



V. DRIVING ACROSS THE FAN 



We now consider the case where the fan of the null is sheared rather than the spine. We 
re-run the simulations with B = B (x, —2y, z), and again drive in the y-direction on the 
x-boundaries. This time we use a uni-directional driving profile, which has the disadvantage 
of compressing the plasma at the boundaries, causing a few extra numerical difficulties, but 
has the advantage of shearing the fan plane in the same direction over the whole 2/z-plane 
for each x-boundary. Specifically, we take 

v = V (t)ir (1 - tanh 2 (A^/H)) (l - tanh 2 (A^/^)) y, (2) 

where Vq is again defined by Eq. ([T]). We take v$ = 0.02, r = 1.8, A y = 12, A z = 5, domain 
dimensions X\ = 0.5, Y\ = Z\ = 3, (3 = 0.05 and Bq = 2 (so that the travel time for the 
disturbance, which propagates at the Alfven speed, to reach the null is approximately the 
same as in the spine shearing cases). 

The evolution of the null point for an ideal monatomic gas (7 = 5/3) is very similar to the 
case where the spine is driven. Once again the disturbance focuses towards the null point, 
this time along its fan, and drives it to collapse. A current sheet which spans the spine and 
fan results [Fig. EJ^a)]. This is expected by comparison with the behaviour of wave-like shear 
perturbations^ 1 ^. However, in the incompressible analytical solution of Craig & Fabling^, 
a shear of the fan leads to tubular current structures aligned to the spine of the null. 
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FIG. 6: (Colour online) Current density |J| in the z = plane, at the time of its temporal peak, 
for (a) 7 = 5/3, and (b) 7 = 10, for the case of driving across the fan. 

Examining the behaviour for larger values of 7, we find that compressibility seems to 
have a similar effect to the spine driving case, but spine current sheets do not develop. 
Specifically, decreasing the compressibility again means that the null does not collapse to 
the same extent, though rather than spreading along the spine as predicted by the analytical 
solutions, the current again spreads along the fan [Fig. 0(b)]. 

These results provide strong evidence that spine current sheets are not dynamically ac- 
cessible, at least in the absence of strong (super-Alfvenic) inflows to drive the localisation. 
This result has previously been anticipated by Titov et a/.—. We instead expect tubular 



spine current structures to be associated with rotational motions, see Ref. [34j. Within 
these tubular structures, we expect the current to flow parallel to the spine, corresponding 
to field lines spiralling around the spine. By contrast, the current in the incompressible 
'spine current' solutions^ is directed parallel to the (undisturbed) fan plane (while being 
localised close to the spine). 



VI. SUMMARY 



We have presented the results of 3D resistive MHD simulations of a driven 3D null point. 
We focussed on the effect of moving from a compressible plasma towards an incompressible 
one, by varying the ratio of specific heats, 7, in our simulations. This was found to strongly 
affect the resulting current sheet formation, both qualitatively and quantitatively. 
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We considered first the case where the spine of the null is sheared from the boundaries. 
For an ideal, monatomic plasma (7 = 5/3, compressible), the spine and fan of the null 
collapse towards one another, and a strongly focused current sheet forms at the null, locally 
spanning the spine and fan. However, as 7 is increased, the collapse of the null, and in 
particular of the fan plane, is suppressed. The current sheet spreads increasingly along the 
fan surface, which remains increasingly planar throughout the simulation runs. In addition, 
rather than forming a stagnation point flow as the null collapses, the plasma flow within 
the domain stays approximately parallel to the planar fan surface for large 7. The same 
effect was found when j3 was increased rather than 7, due to the physically similar nature of 
increasing either parameter, as discussed earlier. Quantitatively, the peak current and peak 
reconnection rate both drop significantly as 7 (or (3) is increased (see also Ref. 4j). 

Considering the case where the boundary shearing was applied across the fan plane of 
the null rather than the spine, we found similar behaviour. In particular, the null point 
collapse is suppressed, and a more spatially diffuse current structure is found, localised to 
the fan surface. Our results provide strong evidence that the steady-state analytical fan 
current sheet solutions of Craig et a/.— are in fact dynamically accessible in a fully 3D, 
incompressible plasma. However, they also imply that the equivalent spine current sheet 
solutions^ are not. Examining the fan current sheet solutions, it appears that the reason 
why they break down in a compressible plasma is the enormous pressure gradients which 
are required to maintain the imposed symmetry. These pressure gradients scale inversely 
with the resistivity, and so in astrophysical plasmas become unphysically large. 
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